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Abstract 

We estimate the electromagnetic effect of the spectator charge on azimuthal anisotropics observed 
in heavy ion collisions. For peripheral Pb+Pb reactions at the top energy of the CERN Super 
Proton Synchrotron, -y/sjvw = 17-3 GeV, we predict this effect to bring very large distortions to 
the observed directed flow, vi, of positive and negative pions emitted close to beam rapidity. The 
overall magnitude of this effect is comparable to values of vi reported by the WA98 experiment. 
We argue that also at lower rapidities, the spectator induced electromagnetic effect may result in 
the splitting of values of vi observed for positive and negative pions. Such a splitting is visible in 
the data reported by the STAR Collaboration from the RHIC Beam Energy Scan. 

Both effects are sensitive to the space-time scenario assumed for pion emission. Therefore, they 
bring new information on the collision dynamics. 
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I. INTRODUCTION 



Azimuthal correlations between particles and the reaction plane, observed in non-central 
collisions, constitute one of the main subjects of heavy ion collision studies. The primary 
nucleon-nucleon collisions being believed to be insensitive to the reaction plane, the latter 
phenomena are considered to be the source of information on collective effects present in 
the heavy ion reaction. A rich phenomenology [l|-[5| has therefore been developed. Among 
others, the success of the hydrodynamical approach in describing these phenomena is to 
be quoted here. 

The azimuthal correlations are usually quantified in terms of the Fourier coefficients of 
the azimuthal distribution of the outgoing particles w.r.t. the reaction plane. These are 
defined by 0]: 

Vn = (cos[n(0 - ^,)]) , (1.1) 

where (f) denotes the azimuthal angle of the emitted particle, while \E'r is the orientation of 
the reaction plane defined by the impact parameter vector b. Specifically, the first-order 
coefficient 



t;i = (cos(0 - ^^)) (1.2) 

refiects the sidewards collective motion of the particles and is known as directed fiow. A 
very sizeable amount of experimental data on directed flow is at present available from the 
low energy regime up to the LHC It is known that Vi depends on collision centrality. 



particle type, transverse momentum and rapidity. For instance, the rapidity dependence of 
pion directed flow at SPS energies is known to consist of a smooth passage through zero 
at mid-rapidity p^], while large values of vi, of the order of 0.2-0.25, are reported in the 
vicinity of target rapidity, in semi-central collisions [13]. 

In the present paper we point at the possible electromagnetic component of directed flow. 
We estimate the size of the spectator-induced electromagnetic interaction on the collective 
motion of positive and negative pions. The spectator-induced electromagnetic effect on 



charged pion ratios (f^) was discussed in our earlier paper 18 1. This effect is in fact 
observed in experimental data [19], and results in very large distortions in ratios of charged 
pions observed at high rapidities and low transverse momenta. Measurements of directed 
flow offer a new way of spotting the influence of the spectator charge on the emission of flnal 
state particles. They also provide a new source of information on the space-time evolution 
of the non-perturbative process of pion production^. 

Our present analysis will be focussed on peripheral Pb+Pb collisions occurring at the top 
energy available to SPS heavy ion experiments {^/sNN = 17.3 GeV). We choose peripheral 
collisions because they are characterized by the largest spectator charge and also large 
values of directed flow. We focus on the SPS energy range by virtue of the abundance 
of experimental data on the rapidity dependence of Vi, which also gives the possibility of 
comparison with our results as explained in Section IIIII We note that the energy range of 
the SPS partially overlaps with that of the RHIC Beam Energy Scan, addressed in short in 
Section |Vl 



^ The cited paper [18| also contains a description of various works made in the past on electromagnetic 



effects in low and high energy nuclear collisions. For more details on this subject, see 20l437|. 



2 



b 




spectator (L) 



FIG. 1: (Color online) Our model of the peripheral Pb+Pb collision. The two spectator systems 
are denoted left (L) and right (R), respectively. 

This paper is organized as follows. Section [TTl contains the description of our Monte Carlo 
model. In Section UTTl we present the results of our simulation, together with the comparison 
to experimental data. The discussion of our results is made in Section [TVl In Section |Vl we 
comment on the experimental data from the RHIC Beam Energy Scan. Our conclusions are 
presented in Section IVll 



II. THE MODEL 



A detailed description of our model can be found in [18[. Only the aspects that are 
relevant for the present analysis will be listed below. Our aim is to obtain a realistic estimate 
of the influence of the spectator- induced electromagnetic interaction on the directed flow of 
pions. On the other hand, we wish to avoid the detailed discussion of the complex and 
poorly known mechanism of soft particle production. Therefore we decide on a maximally 
simplified approach: 

(a) We assume, as an example, a peripheral Pb+Pb collision involving Upart = 60 partic- 
ipating nucleons at ^/sNN = 17.3 GeV. This is shown in Fig. [H The two spectator 
systems are modelled as two homogeneous, Lorentz-contracted spheres. The reac- 
tion plane is defined by the collision axis (dashed line in the figure) and the impact 
parameter vector b. 

(b) Charged pions are assumed to be emitted from a single point in position space, that is, 
the original interaction point. The time of pion emission t^; is a free parameter of our 
model; this parameter sets the initial conditions for the electromagnetic interaction. 
For peripheral Pb+Pb collisions studied here, the initial spectra of charged pions are 
assumed to be similar to those in nucleon-nucleon collisions. 

(c) Charged pions are then numerically traced in the electromagnetic field of the spectator 
charges until they reach the distance of 10,000 fm away from the interaction point and 
from each of the two spectator systems. The fragmentation of the spectator systems 
is neglected. We do not consider the effects of participant charge, strong final state 
interactions, etc. 
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Several clarifications should be added to the above. 



The simplified geometry of the peripheral collision presented in Fig. [T] was determined 
on the basis of a dedicated study, discussed in detail in [ll]. A geometrical Monte 



Carlo simulation served this purpose. This used proper nuclear density profiles [38 
and assumed the elementary nucleon-nucleon cross-section equal to 31.4 mb, in good 
agreement with experimental p+p data at this collision energy [39[. The geomet- 
rical impact parameter corresponding to 60 participating nucleons was found to be 
bgeom = 10.61 fm. Additionally, the center of gravity of each of the two spectator 
systems was found to be displaced by Ab = 0.76 fm relative to that of the original 
nucleus. The average spectator charge was Q = 70 elementary units. 
Considering the exact shape of the spectator system as unimportant for our subse- 
quent studies, and for the sake of clarity, we modelled the two spectator systems as 
homogeneously charged spheres. In the rest frame of each sphere, its density was the 
standard nuclear density p = 0.17/fm^. The center of each sphere was additionally 
displaced by 0.76 fm in order to match the center of gravity of the spectator system. 
Thus our effective impact parameter (distance of closest approach between the two 
spheres' centers) was b = bgeom + 2A6 = 12.13 fm. 

(b) The simplification of initial conditions for pion emission (reduction of the emission 
zone to a single point in space and time, as discussed above) gives a convenient way 
to estimate the sensitivity of the electromagnetic effect to the basic characteristics of 
pion production (the pion formation time, the distance between the pion formation 
zone and the two spectator systems, etc). The initial momentum spectra of pions 
(before the action of the electromagnetic field) are assumed to obey wounded nucleon 
scaling and to be similar to underlying nucleon-nucleon events. As such, they are 
described by an analytical parametrization of average pion ( ^^~^^ ) distributions in 



p-|-p collisions, recorded by the NA49 experiment at the SPS |39|. More details on this 
parametrization are given in [isj. For simplicity, distributions of positive and negative 
pions are assumed to be identical. The experimental data cited above are expressed 
in terms of the Feynman variable^ Xp = 2pL/y/s and of transverse momentum pt, 
and cover the region from xp = to 0.85 and from p^- = to 2.1 GeV/c. Only a 
small extrapolation towards higher xp needs to be applied in the present analysis. 
The uncertainty of this extrapolation has little or no effect on the results presented in 
Section lllli It should be underlined that full azimuthal symmetry is assumed for the 
initial emission of pions. 

Charged pions, with their initial momentum vector defined above in point (b), and 
weighted by their initial distribution, are subjected to the electromagnetic field of the 
two spectator systems. The pion trajectory fv is given by the classical relativistic 
equation of motion: 

^ = vAr,t)= (2.1) 



All the kinematical variables addressed in this paper will be defined in the cm. system of the collision. 
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where m^r is the pion mass, and the pion momentum is defined by the Lorentz force 
acting on the pion: 

^ = F^{f, t) = (E{r, t) + v^{r, t) x 5(f, t)) . (2.2) 

Here, q-,; is the pion charge, while E{f,t) and B{f,t) are standard superpositions 
of fields from the two spectator systems. Our equations take account of relativistic 
effects, including retardation. Technically, the propagation of the pion is performed 
numerically by means of an iterative procedure made in small steps in time, with 
variable step size. The iteration proceeds until the pion is 10,000 fm away from the 
interaction point and from each of the two spectator systems in their respective rest 
frames. Negative pions which do not escape from the potential well induced by the 
spectator system are rejected and do not enter into our final state distributions. 



III. RESULTS 

In this Section we present the results of our analysis. These will be displayed as a function 
of the scaled rapidity y/yheam, where y and yheam are the rapidity of the pion and of the 
incoming nucleus in the collision cm. system. We decide on this particular variable in order 
to simplify the interpretation of our results, and for an easier possible comparison with other 



collision energies (see also the discussion made in 12|, |41| ) 



A. Electromagnetic effect from one and two spectators 

We start by discussing the basic features of the spectator-induced electromagnetic effect 
on the directed fiow of pions. In particular, this will include the role of each of the two 
spectator systems in the overall effect. Fig. [2] shows the directed flow vi of vr"^ and vr", 
induced by their electromagnetic interaction with the spectator charges. The value of Vi is 
integrated over the transverse momentum of the pion, from = to 1 GeV/c. All the 
simulations presented in the flgure assume the simplest situation where the pion emission 
time t^; is equal to zero (immediate pion creation). 

It should be underlined that as our simulation contains no initial azimuthal anisotropy 
(Sect, ini), any non-zero value of Vi apparent in Fig. [2] is indeed solely due to the elec- 
tromagnetic interaction with the spectator systems. Therefore, we will refer to it as 
"electromagnetically-induced directed flow" in order to differentiate it from the "standard" 
flow phenomena which are caused by the strong interaction. In this context, several remarks 
are in order. 

1. As it is clearly apparent from the flgure, the electromagnetic fleld induced by the 
spectator charges exerts a strong influence on the directed flow of charged pions. It 
can be held responsible for large values of vi, exceeding 0.2 (20%) for negative pions 
close to beam rapidity. 

2. The directed flow induced electromagnetically by the two spectators (red solid curve, 
which goes through the black dots in the flgure) exhibits a characteristic structure as a 
function of rapidity. Qualitatively, this structure recalls the overall behaviour of pion 
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FIG. 2: (Color online) Spectator-induced electromagnetic effect on directed ffow of (a) positive and 
(b) negative pions, in peripheral Pb-|-Pb collisions at ^/snn = 17.3 GeV. The green solid curve 
shows the directed flow induced electromagnetically by the right (R) spectator. The blue solid 
curve shows the directed flow induced electromagnetically by the left (L) spectator. Black dots 
show the result of the addition of these two curves. The red solid curve displays the result of the 
simulation including both spectators. Note: all the simulations assume Ie = fm/c; the blue solid 
curve is obtained by reflection of the green solid curve about midrapidity. 



directed flow as it is reported from experimental data (see e.g. (l2|, ll7|). The directed 
flow displays a smooth transition through zero at y/yheam = 0, with large values of 
flow (anti-flow) at large negative (positive) rapidity. 

3. The sign of electromagnetically-induced directed flow for negative pions is opposite 
with respect to positive pions. For absolute values of vi, differences in the shape of 
the two curves are seen close to beam rapidity. 

4. What specifically follows from the above is that the spectator-induced electromagnetic 
interaction may result in the splitting of values of Vi observed for positive and negative 
pions. This splitting exhibits a strong rapidity-dependence, and reach very large values 
close to beam rapidity. 

5. The directed flow which is induced electromagnetically by each of the two spectators 
separately is also presented in the figure (the two spectator systems are denoted "left" 
and "right" as in Fig. [T]). The contribution of a single spectator does not remain 
confined to its hemisphere of the collision. On the contrary, it extends into the opposite 
hemisphere; thus the "right" spectator moving at positive rapidity exerts its influence 
on the pion vi down to and beyond y/yheam = —1, and vice-versa for the "left" 
spectator. 

6. The value of vi obtained as the result of the electromagnetic interaction with both 
spectators (red solid curve, which goes through black dots in the figure) appears equal 
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FIG. 3: (Color online) Spectator-induced electromagnetic effect on directed flow of positive and 
negative pions, in peripheral Pb+Pb collisions at ^/snn = 17.3 GeV. The top and bottom panels 
correspond to different values of the pion emission time tg. 



to the result of direct addition of the two single-spectator curves (displayed by the 
black dots). 

Thus, as it it becomes evident from the above, the spectator-induced electromagnetic 
interaction appears strong enough to result in sizeable effects, and to play an important 
role in flow phenomena involving charged pions in heavy ion collisions. This imposes the 
necessity of more detailed studies which will be presented below. 



B. Dependence on initial conditions 

The central issue of this paper is the sensitivity of the spectator-induced electromagnetic 
effect to the space-time evolution of the collision. This will be addressed in the present 
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Section. 

Fig. [3] shows the rapidity-dependence of the electromagnetically-induced directed flow of 
positive and negative pions, obtained assuming two different values of the pion emission 
time: = and 1 fm/c. The values of Vi are obtained by integration over px of the pion, 
with integration limits defined as in Section IIII Al above. The curves in panels (a) and (b) 
of the figure are the same as the respective curves in Fig. [2J The following remarks are to 
be made. 

1. The electromagnetically-induced directed flow displays a clear dependence on the pion 
emission time. In the region of beam rapidity, the position of the minimum (maximum) 
in the valley (peak) in pion vi shifts from y/uheaxa = 1 at t^; = fm/c to above 
y/yhea.m = 1.2att£; = lfm/c. At the collision energy discussed here, this shift in 
y/Vheaui corresponds to a shift of about 0.6 units in rapidity. The actual shape of the 
valley (peak) also exhibits sensitivity to Ie- 

2. Also in the region of lower rapidities, and down to midrapidity, the change mtE results 
in changes of Vi. For instance, at y/ybcam = 0.5 the absolute values of Vi go from about 
2.4% at t^; = fm/c down to about 0.5% at t^; = 1 fm/c. 

3. What follows from the above is that the electromagnetically-induced splitting of vi for 
77+ and 7r~, mentioned in Section [III At will also exhibit sensitivity to Ie- 

These observations demonstrate that the spectator-induced electromagnetic effect is in- 
deed sensitive to the initial conditions, set for the electromagnetic interaction by the emission 
time Ie- Changes in the pion emission time, or equivalently in the position of the pion for- 
mation zone with respect to the two spectator systems, will result in changes of vi which 
should be observable in experiment. In other terms, the electromagnetic effect on directed 
flow is indeed sensitive to the space-time evolution of pion production, and provides new 
information on the dynamics of the collision. 

C. Dependence on transverse momentum 

Up to now, for simplicity, only the rapidity-dependence of pT-integrated directed flow was 
discussed in this paper. In the present Section we study its dependence on the transverse 
momentum of the pion. This is shown in Fig. HI The electromagnetically-induced directed 
flow exhibits a strong dependence on pion pt- In particular, this is valid close to beam 
rapidity where very large absolute values of Vi (of the order of 0.6 for positive pions) are 
attained at the lowest considered value of transverse momentum, px = 75 MeV/c. The 
absolute value of the electromagnetically-induced vi rapidly decreases with increasing pt- 
On the other hand, little or no dependence on pion px is apparent in the region of lower 
absolute values of y/yheam, down to midrapidity. The overall pattern of dependence on initial 
conditions, discussed in Section UlI Bl above, repeats itself also at fixed values of pt- 

D. Comparison to experimental data 

Taking into account the charge-asymmetric nature of the electromagnetic effect discussed 
here, any comparison of our work to experimental data must involve measurements of di- 
rected flow of pions of a given charge, either tt"*" or n". This constitutes an important 
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FIG. 4: Spectator-induced electromagnetic effect on directed flow of positive and negative pions 
in peripheral Pb+Pb collisions at y^ijvw = 17.3 GeV, shown at fixed values of pion transverse 
momentum: pr = 75 MeV/c (solid), = 125 MeV/c (dashed), pr = 175 MeV/c (dotted). The 
top and bottom panels correspond to different values of the pion emission time tE- 



practical difficulty as most of results published on the rapidity-dependence of f i involve (if 
particle identification is available at all) only "charged pions", i.e., summed vr"*" and vr". One 
data set which fulfills our requirement comes from the WA98 experiment, and contains a 
measurement of positive pions near target rapidity made by means of the Plastic Ball detec- 
tor, in Pb+Pb collisions at ^/sNN = 17.3 GeV. This measurement is presented in but a 
more detailed description of the analysis can be found in j42[. The latter reference specifies 
the corresponding centrality definition as 40-80% of the Pb-|-Pb cross-section. Following the 
discussion made therein, we expect that differences w.r.t. our definition of the peripheral 
collision (Section |Tl]) would only have a small effect on the actual values of fi. 

Fig. [5] shows the experimental data from WA98, superimposed with the rapidity- 
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FIG. 5: Comparison between experimental data on directed flow of positive pions obtained 



by the WA98 experiment (the data points are redrawn from 17|]), and our simulation of 
electromagnetically- induced directed flow of vr^. The three curves correspond to three differ- 
ent values of the assumed pion emission time tE- Note: the data points at positive y/?/beam are 
obtained by reflection about midrapidity. 



dependence of the electromagnetically-induced directed flow as obtained from our work. 
For the latter, three values of the pion emission time are assumed: tE = 0, 0.5, and 1 fm/c. 
Our results are integrated over px from to 1 GeV/c. The curves corresponding to t^; = 
and 1 fm/c are the same as in Fig. El 

As it can be seen in the figure, the values of Vi obtained from our simulation, and resulting 
exclusively from the spectator-induced electromagnetic interaction, appear comparable to 
those measured by experiment. 

The above statement should clearly be taken with some caution. Account should be taken 
of the simplicity of our model, and of various detailed issues related to the measurement 17 . 
These problems are beyond the scope of the present paper. 

Nevertheless, it seems clear that a very sizeable part of the measured directed flow of 
7r~^ may come from the electromagnetic origin. We note that this observation is consistent 
with our analysis of the influence of spectator charge on tt^/tt^ ratios in peripheral Pb+Pb 
collisions at the same energy. Here, pion emission times assumed in the range = 0.5- 
1 fm/c also result in good agreement between our model and experimental data |43|. 



IV. SUMMARY AND DISCUSSION 

Summing up the principal results of our work, the following can be said: 
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1. We devised a simple model of the Pb+Pb collision (Section [IT])- This model does not 
contain any initial azimuthal anisotropy in pion emission but includes the electromag- 
netic interaction between the pion charge and the two spectator systems. 

2. This model predicts sizeable values of vi (electromagnetically- induced directed flow) 
for positive and negative pions in the final state of the collision. For positive pions, 
the predicted values are comparable to those seen in experimental data. As such, 
a very large part of positive pion directed flow may come from the electromagnetic 
interaction. 

3. The sign of this effect is opposite for positively and negatively charged pions. This 
will result in a splitting of values of vi observed for vr^ and vr". 

4. The size of electromagnetically-induced directed flow, and therefore also the size of the 
electromagnetic splitting addressed above, depends on the pion emission time ts- The 
latter parameter defines, in our model, the details of the space-time evolution of pion 
production, like e.g. the time of pion formation, or the position of the pion formation 
zone with respect to the spectator system(s). 

5. As it follows from the above, the electromagnetically-induced directed flow, and the 
resulting splitting of Vi for tt"*" and vr", bring new information on the space-time dy- 
namics of the Pb-hPb reaction. 

This has, in our view, important consequences for future studies of directed flow. It is 
of course clear that our model contains important simplifications. Among others, it inten- 
tionally neglects the initial pion flow induced by the strong rather than the electromagnetic 
force^. The inclusion of such phenomena is the domain of more detailed models and we leave 
it for future studies. However, it is clear from our work that the electromagnetic splitting 
of fi(7r"'") and fi(7r~) addressed above will remain present also in more sophisticated phe- 
nomenological descriptions of the heavy ion collision. This inspires a set of proposals and 
remarks listed below. 

• The study of rapidity-dependence of directed flow made separately for pions of a given 
charge, i.e., vr"*" and vr", will bring additional information on the dynamics of the 
heavy ion collision. This information will add up to the already rich phenomenological 
content of summed charged pion flow studies. The electromagnetic splitting of vi 
for positive and negative pions will provide information about the position of the pion 
formation zone with respect to the spectator systems. Specifically, also the lack of such 
splitting at a given rapidity would imply that the formation of pions occurs relatively 
far from the nearest spectator. 

• The situation should be similar for other produced particles, like e.g., charged kaons. 
The situation for final state protons should be at least partially different as protons at 



forward rapidity will come also from spectator fragmentation [20|, 14^, |45|. However, 
also in this case the role played by the electromagnetic repulsion of protons from the 
spectator systems should be investigated in detail. 
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We keep in mind that the inclusion of the latter is necessary in order to describe the summed charged 
pion data (see, e.g., [l^). 
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• Last but not least, our analysis brings attention to asymmetric nuclear collisions, be it 
A+B, d+A, or p+A reactions. In such collisions, the electromagnetic influence of the 
two spectators would not cancel out at midrapidity, and therefore the role of single 
spectators could be isolated, yielding further information on the space-time evolution 
of the collision (see Fig. [2] for comparison). 

All of the above observations are consistent with the results of our study of spectator- 
induced electromagnetic forces acting on charged pion ratios [18]. However, we feel that the 
precision of the present experimental measurements of the size of directed flow, the key role 
played by flow measurements in the present phenomenology of heavy ion collisions, and the 
overall interest of the community in such studies, make the consideration of the role played 
by electromagnetic interactions important for future analyses. 



V. DATA FROM THE RHIC BEAM ENERGY SCAN 

In the context of the discussion made above, and based on the comparison between the 
results of our simulation and existing SPS experimental data, it is important to note that 
the RHIC Beam Energy Scan (BES) program can also address the search of the spectator- 
induced electromagnetic effects discussed here. 

Specifically, preliminary data on directed flow of positive and negative pions from the 



STAR experiment have already been reported in [46|. The measured values of Vi{7r~^) 
and Vi{7r~) indeed display a signal of splitting which increases with increasing rapidity. 
At positive rapidity, Vi{tt~^) remains below vi{tt~), which follows our predictions for the 
electromagnetically-induced directed flow, formulated in Sections IIII Al and IIII Bl above. 

Presenting final conclusions on the nature of this effect probably requires more precise 
modelling of the space-time mechanism of charged pion production, including also its cen- 
trality dependence. However, we note that the overall magnitude of the splitting of vi values 
seen by experiment looks (roughly) comparable to that predicted in the preceding Sections, 
for the largest assumed values of the pion emission time tE and in the relatively limited 
range of y/j/beam available to the measurement j46|. Thus, we consider that the observed 
splitting is indeed a signature of electromagnetically-induced directed flow. 



VI. CONCLUSIONS 

The electromagnetic field induced by the presence of spectator systems may result in 
sizeable distortions of the directed flow pattern observed for pions of a given charge (that is, 
TT"*" and TT") produced in heavy ion collisions. This implies the presence of the electromagnetic 
splitting of vi for particles of opposite charges. The size of this splitting depends on the 
space-time scenario of pion production and therefore brings new information on the collision 
dynamics. The comparison to experimental data suggests that a very large part of directed 
flow observed for positive pions may come from the electromagnetic origin. 

New experimental and phenomenological studies are needed in order to elucidate these 
questions, and possibly bring new insight into the evolution of the heavy ion reaction 
in space and time. On the experimental side, this implies measurements of particles of 
specific charges, for symmetric and asymmetric nuclear reactions. On the theoretical side, 
the possibility of verification of existing models of the heavy ion reaction by taking into 
account spectator-induced electromagnetic effects should be investigated. Furthermore, 
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studies of the role of spectator-induced electromagnetic interactions for higher harmonics 
of the azimuthal distribution would also be highly indicated. 
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